""" Программа для вычисления минимального временного шага
    расчета для явной динамики в методе КЭ 
    для материала с модулем упругости E и плотностью rho
    с минимальным ребром конечного элемента lengthOfMinFE
    
    E - модуль упругости материала, Па;
    rho - плотность материала, кг/м3;
    lengthOfMinFE - минимальная длина КЭ в сетке;
    C - скорость звука в материале, м/с;
    timeOfMinFE - время прохождение упругой волны через минимльную длину КЭ в сетке.
    """

import math
E = 2e11
rho = 7850

C = math.sqrt(E/rho)
print(f"Elastic constant: {E} \nDensity: {rho} \nSpeed of sound: {C}")

lengthOfMinFE = 1
timeOfMinFE = lengthOfMinFE / C
print(f"Minimum timestep for calculation: {timeOfMinFE}")